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Abstract 

The next-to-next-to-leading order (NNLO) non-singlet QCD analysis of the exper- 
imental data of the CCFR collaboration of Ref.[l] for the xF 3 and F 2 structure 
functions of the deep-inelastic scattering of neutrinos and antineutrinos on the nu- 
cleon by means of the Jacobi polynomial expansion method is made. The target 
mass corrections are also taken into account. We demonstrate that the NNLO 
corrections have the important effect: they are decreasing the difference between 
the values of a s (Mz), extracted from the NLO fit of xF 3 data and from the NLO 
non-singlet fit of the combined xF$ and F 2 CCFR data. The obtained NNLO re- 
sults for a s (M z ) are : a s (M z ) = 0.109 ± QM3(stat) ± 0.005(syst) ± 0.003(theor) 
(from xF 3 data) and a s (M z ) = 0.111 ± 0.002(stat) ± 0.003(syst) ± 0.003(theor) 
(from the combined xF 3 and F 2 data). We also estimate the values of the order 
0(a 2 )-corrections to the Gottfried sum rule. 



1. The recent measurements of the CCFR collaboration provide the most precise up to now 
experimental results for the structure functions (SFs) xF s (x,Q 2 ) and F 2 (x,Q 2 ) of the deep- 
inelastic scattering (DIS) of neutrinos and antineutrinos on nucleons [|J . The next-to- leading 
order (NLO) QCD analysis of the available at present data, obtained at Fermilab Tevatron, 
was made in Ref.|]J with the help of the approach, based on the solution of the correspond- 
ing integro-differential equation 0, ||. The results for the parameter A^, were extracted 
separately from the fits of the observed behavior of the xF 3 SF and from the nonsinglet (NS) 
analysis of the combined available data for xF% SF at x < 0.5 and for F 2 SF at x > 0.5 in 
the energy regions where the higher-twist corrections were considered as the negligible small 
effects. In the process of the latter fits the statement was made that the substitution of the 
data for xF 3 by the data for F 2 at x > 0.5 is improving the experimental precision of the 
results obtained due to higher statistics of F 2 in this region. 

The value of the parameter A^-, which resulted from the NLO fits of the data for xF 3 SF, 

A4 o 



is n 



Ag) = 171 ± 32(stat) ± 5A(syst) MeV 



MS 

while the combined NLO NS analysis of xF 3 and F 2 SFs data gave [I 

A^L = 210 ± 28{stat) ± Al(syst) MeV . (2) 

Let us remind that there is also another method of the analysis of DIS data, which is based on 
the idea of the reconstruction of the SF from its Mellin moments using the expansion of the 
SF into the series over the Jacobi polynomials 0]. This method was developed in Ref. 0] and 
Refs.|6|, [jj and successfully applied in the process of the NLO fits of DIS data of the BCDMS 



collaboration || and the mentioned above data of the CCFR collaboration P, |I0|. 

In the last several years the considerable progress was achieved in the area of the analyt- 
ical evaluation of the characteristics of DIS at the next-to-next-to-leading order (NNLO) of 



perturbative QCD in the MS'-scheme. Indeed, the NNLO corrections were found in Ref.JTT| 
for the coefficient functions of the SF F 2 , while in Refs.[^, [13|| the similar NNLO calculations 
were made for the SF xF%. The result of Ref.fT2|| is in agreement with the one for the NLO 



corrections to the n = 1 moment of xF%, namely to the Gross-Llewellyn Smith sum rule |14 



while the results of Ref. [11] are in agreement with the expressions for the NNLO corrections 
to the coefficient functions of the NS moments of F 2 SF, calculated analytically in the case 



of the even moments with n = 2,4,6,8,10 [Jl~5 |. Moreover, in the case of n = 2,4,6,8 the 



analytical expression of the NNLO corrections to the anomalous dimensions of the NS mo- 
ments of F 2 SF are known at present flEfl . Quite recently the analytical result for the NNLO 
correction to the anomalous dimension of even n = 10 NS moment of F 2 SF became available 
[I7|. Together with the expression for the NNLO correction to the QCD /5-function |I8|, which 



was confirmed in Ref.|T9[, the results of Refs.fl l||-fl7|j are forming the theoretical background 



for the extraction of the value of the parameter A^L from the NS fits of the DIS data at the 
level of new theoretical precision, namely with taking into account the effects of the NNLO 
corrections. 

We think that it is still rather difficult to push this program ahead in the framework of the 
integro-differential formalism. Indeed, the analytical dependence of the NNLO corrections to 
the anomalous dimensions from the number of the related moments n is unknown at present. 
Moreover, in view of the definite limitations of the possibilities of the distinguished methods 
of Ref. f2Cfl , used in the calculations of Refs.[[l4|-[p!7|, this information can not be 

obtained without additional theoretical efforts []. This fact does not allow one to calculate at 

1 The calculating method, analogous to the ones developed in the first two works of Ref. (2(j] was also 
developed for other purposes in Ref.[f22f 
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present explicitly the NNLO correction to the kernel of the corresponding integro-differential 
equation. 

It should be stressed that in the framework of the method of the Jacobi polynomials it is 
necessary to use the information about the renormalizat ion-group evolution only of the fixed 
number of moments. To study the evolution of the NS moments at the NNLO level it is possible 



to use the results of the exact calculations of Refs.(TT]]-[|T3j and to supplement the information 
about the NNLO coefficients of the NS anomalous dimensions of the certain even moments 
[ IBf by some smooth interpolation procedure [23]. In Ref. [23 this program was realized for 



the extraction of the value of the parameter Ajg^ from the NS fits of the data for the F 2 SF, 
provided by the BCDMS collaboration, in the approximation when the effects of the target 
mass corrections (TMC) were neglected. Note, however, that this analysis is non-complete, 
since in the case of the NNLO fits of the data for the DIS of charged leptons on nucleons it 
is also necessary to take into account the information about the NNLO QCD corrections to 
the singlet coefficient functions and singlet anomalous dimensions, recently calculated in Ref. 

This paper is devoted to the NNLO NS QCD analysis of the available data provided by 
the CCFR collaboration [|TJ. First, in order to compare the outcomes of the applications 
of the Jacobi polynomials expansion method with the NLO results, obtained by the CCFR 
collaboration, we are performing the detailed NLO fits with taking into account the TMC 0. 
Then we are performing the detailed NNLO fits of the available CCFR data. Using the certain 
interpolation/extrapolation procedure we are also presenting the estimates of the previously 
unknown order 0(a;?)-corrections to the Gottfried sum rule. 

2. We first discuss the theoretical basis of the method used in this work. Let us define 
the Mellin moments for the NS SF xF 3 (x,Q 2 ): M^ S (Q 2 ) = $ x n-1 F 3 (x, Q 2 )dx where n = 
2,3,4,.... The theoretical expression for these moments obey the following renormalization 

group equation (pJL + f3(A s )^ - 7 JS(A a )) M» s \Q 2 / )j? , A s (fi 2 )) = where A s = a s /(4n). 

For the simplicity we consider the case when the factorization scale and the scale of the 
ultraviolet renormalizations are equal to each other. The renormalization group functions are 
defined as ^ = f3(A s ) = -2£ 4 > A^ +2 , /x^gf^ = >y%l(A.) = E^oT&W^ 1 where 
Z^ s are the renormalization constants of the corresponding NS operators. The solution of 
the renormalization-group equation can be presented in the following form 



M^(Q 2 ) 
M» S (Q 2 ) 



exp 



MQ>) lP s (x) dx \ ctil(A s (Q 2 )) 
MQ 2 ) P( x ) 



Cns(A s {Qo)) 



(3) 



where M^ s (Ql) is the phenomenological quantity related to the factorization scale dependent 
factor. It can be parametrized through the parton distributions at fixed momentum transfer 
Qo- I n our studies we will consider one of the possible forms of the parametrization of xF$ 
SF, choosing M* S (Q 2 ) = ^ x n ~ 2 A{Q 2 )x b ^ {1 - x)< Q l\l + -f(Q 2 )x)dx with 7 ^ 0. 

At the NNLO the expression for the coefficient function has the following general form 
C$l(A a ) = 1 + C (1 \n)A s + C^(n)Al while the corresponding expansion of the anomalous 
dimensions term in Eq.(|3|) is 



exp 



(A s (Q 2 ))^y^ [1+p{n)As{Q 2 ) + g(n)(As( g2 ))2] (4) 



In the similar analysis of Ref. M [lOfl the effects of the TMC were neglected. 



where p(n) and q(n) are denned as p(ri) = \ Ims^H _ 7 ^|( n ) & _ I[2p(n) 2 + 7 ^| w -f 



/3i Pa 



— /fi^ lNsi 71 ) ~~ Tats ( n ) /I?] • Notice that in our normalizations /? , Pi and [3 2 read /3 — 
11 - (2/3)/, (3i = 102 - (38/3)/, (3 2 = (2857/2) - (5033/18)/ + (325/54)/ 2 . 

Now we are ready to present the basic formula for the Q 2 -evolution of the NS-moments, 
which will be used in our fits: 

M^ S (Q 2 ) _ ^ s (g 2 )l^H/ 2 ^[l + p(n)A s (Q 2 ) + gH(^(g 2 )) 2 ]C'S(^(g 2 )) 



A S (Q 



[l+p(n)A s (Q 2 ) + q(n)(MQl)) 2 ] C^ S (A S (Q 2 )) 



Eq.© can be expressed in terms of the inverse powers of the \n(Q 2 / Aj^) using the tradition- 
ally considered expansion 



' ' U ; 4vr PoHQVAy) f3lln 2 {Q 2 / A 2 ^) 
(3 2 {ln 2 MQVAf^) - lnln(Q 2 /Ay) + ^ - 1) 



+ 



/^(QVA 2 ^ 



The similar analysis can be also done for the NS part of F 2 SF. This analysis can be applied 
in the kinematical region x > 0.5 where the valence quark distributions are dominating over 
the distributions of the sea quarks and gluons. Therefore in this kinematical region the singlet 
nature of F 2 SF can be neglected and the approximation xF 3 = F 2 , considered by the CCFR 
collaboration |TJ], can be used after taking into account the effects of the TMC. It should 
be noted, however, that this procedure has definite theoretical uncertainties, since the NS 
coefficient functions and the anomalous dimension terms of the xF 3 and F 2 SFs are rigorously 
speaking slightly different (see discussions below). 



In order to include the TMC in our fits we introduce the Nachtmann moments p4| of the 
SFs F 2 and F 3 : 

M^(Q 2 ) - Cdx^—F^x g 2 ) (3 + 3(n + l)y + n(n + 2)y 2 ) 

M™?(Q 2 ) = /^^3(x,g 2 ) (1+ ( f + + 2 | )T/) , (8) 

where £ = 2x/(l + V), V = \Jl + AM 2 ucl x 2 / Q 2 and M nud is the mass of a nucleon. These 
expressions could be expanded into a series in powers of M 2 ucl /Q 2 |2~5| . Taking into account 
the order 0(M^ ud /Q A ) corrections, we get 

M™ c m = M W ^(Q 2 ) + ~ 2 1} M2 Q 2 d - MZ(Q 2 ) 0) 
, (n + 2)(n + l)n(n-l)M^ „ s . n , , 



2(n + A)(n + 3) Q 4 n+4V ^ ' Q 



(10) 



Ul + 2 )(n + l)nM* ucL n/rNS 2 nf Kvd s 



We have checked that the the influence of the order 0{ fjr 1 ) terms in Eqs.(9),(10) to the 



Q- 1 



outcomes of the concrete LO and NLO fits is very small. Therefore, in what follows we will 
use only the first two terms in the r.h.s. of Eqs.(9),(10). 

The SF is reconstructed from the corresponding Mellin moments using the following equa- 



tion 



x pN max [ 



x,Q 2 ) 



x 



X 



max Tl 



(11) 



where 0" ,/3 are the Jacobi polynomials and a, j3 are their parameters, fixed by the condition of 
the requirement of the minimization of the error of the reconstruction of the SF (see Refs. || [7j ) . 
In our analysis we are considering the region 6 < N max < 10, which corresponds to the number 
of moments 2 < n < N max + 2, or to 2 < n < N max + 4 if the TMC are included up to order 
0(M^ ud /Q 2 )-terms in the r.h.s. of Eqs.(9),(10). Thus we will be interested in the study of 
the behavior of the moments with 2 < n < 14. 

3. Let us now summarize the available perturbative QCD information for the coefficient 
functions of the NS Mellin moments in the MS'-scheme. The corresponding NLO coefficients of 
the Mellin moments of the NS SF F 2 was calculated in Ref. |[26||. Its exact analytical expression 



reads Cp} (n) 



2[S 1 {n)} 2 +3S 1 (n)-2S 2 (n)- 



2Sx(n) , 3 I 4 |_ _2__q 

n(n+l) ~ n ~ n+1 ~ n 2 



where S k (n) = EL 



The similar term in the expression for the coefficient function of the Mellin moments of the 



n 



c 



F 2 ' 



n 



8 2ra+l 
3 n(n+l) 



2q| . Since 



SF xF 3 receives the following additional contribution C Fs 

we will be interested in the numerical values of the corresponding coefficients of the coefficient 
function C^g{A s ) for the fixed number of the moments, we present them in the numerical form 



(2) 

in Table 1, choosing the number of flavors / = 4. The results for the NNLO coefficients Cp 2 (n) 

(2) 

and C F . £ (n) are obtained from the results of Refs. [11,12] and Ref. [13] respectively taking the 
corresponding Mellin moments P[ 

1 1 ' ' 1 ' " ■ / ■ 1 ; . ; i ,,\ j s decreasing with the increase of 



Notice, that the factor A 



n 



the number of moment. The similar feature also holds at the NNLO level for the difference 



A 



(2) 
C 1 



n 



C!$(n)-C%(n) (see Table 1). 



.(2) 



We will discuss now the contributions to the NS anomalous dimensions of the Mellin 
moments of F 2 and xF% SFs. The expressions for the leading order coefficients r Y^sF 2 ( n ) an< ^ 



„(o) 



lNS,F 3 { n ) coincide with each other and read 7^s(^) = | 

coefficients of the NS anomalous dimensions were first obtained in Ref.p7 
were presented later on in the more simplified form [E? 



4Si(n) 



n(n+l) 



The NLO 
These results 



The NS results of Refs. 137], were 



confirmed by the independent calculation in Ref. |29j . It should be stressed, that starting from 
the NLO the analytical expressions of anomalous dimensions of the NS moments of F 2 (xF 3 ) 
SFs can be obtained from the direct calculations of the corresponding Feynman diagrams 
in the momentum space only in the concrete cases of the even (odd) number of the Mellin 
moments. To find out the values of these anomalous dimensions at any n the procedure of the 
analytical continuation should be performed for the cases of odd (even) number of moments 
0. This procedure can be realized using the results of the considerations of the first work of 
Ref. [21 1 and of Ref.0 (the related analysis was described in more detail in Ref. p0| ) . In the 
cases of F 2 (xF 3 ) SFs the procedure results in the application of the "+" ("— ") prescriptions of 



to 



3 The term proportional to CaCf^(1 — z) in Eq.(15) of Ref.[12] contains the misprint. It should be changed 
— (50/3 — z!42/3). After this change the results for the related moments coincide with the results given in 



Table 2 of Ref. [13]. 

4 In principle, the corresponding calculations in the x-space lead to the complete expressions for the anoma- 
lous dimensions without application of the procedure of the analytical continuation. 
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Ref. 0, which shift slightly up (down) the values of the corresponding anomalous dimensions 
of the odd (even) moments. 

We have noticed, as the authors of Ref. 0, that even though the applications of the pre- 
scriptions mentioned above produce a very small change of the values of the NLO corrections 
to the NS anomalous dimensions, they are extremely important for the reconstruction of the 
SFs using the Jacobi polynomial method. In fact when high moments with n > 10 are involved 
in the polynomial reconstruction, the existence of very large factorial terms in the Jacobi poly- 
nomial expansion spoil the convergence of the fitting procedure if the effects of the analytical 
continuation are not taken into account. 



To save the space we will not present here the final analytical expressions for 7^ F2 



11 



and 77V5_f 3 (^) through the poligamma-functions and some other associated to them functions, 
which are continuous in its argument n. However, in Table 2 we present their numerical 
values for 1 < n < 14, fixing the number of flavors / = 4. We emphasize, that using the 
1/n-expansion (which turns out to be the good approximation starting already from n = 4) 



it is possible to show, that the difference between the final analytical expressions of 7^5 Fi 



n 



and 7w 3 ' 



n 



is rather small, namely A„ 



n 



lNS,F 3 ( n ) 



31]. For n = 2, 4, 6, 8 the concrete numerical values of A,. 



7w» ~ 128/(3n 6 ) +0[l/n 7 
are listed in Table 3. For larger 



in 



even values of n the non-rounded off values of A 7 are less than 10 . 

Assuming that the feature, similar to discussed above, will remain true at the NNLO also, 
we will use the calculated in Ref . [|16| , |17| values of j NS F 2 {n) in the NNLO analysis of xF% data. 



In analogy with the considerations of Ref.[[23[], the numerical values of Tns F 2 ( n ) 
for n = 3, 5, 7, 9 has been obtained from the results of calculations of Ref. fll6|, p~7f| (made for 
the case of n = 2, 4, 6, 8 in Ref.fTBfl and of n = 10 in Ref.|i7|) using the smooth interpolation. 



The concrete numerical expressions for the coefficients 7jvs( n ) are presented in Table 2 for 
2 < n < 10. It should be stressed, that in order to verify the validity of the interpolation 
procedure used by us, it is of real interest to calculate the values of 7^ Fs (n) analytically in 
the case of n — 3, 5, 7, 9. Moreover, it is also rather desirable to learn the explicit expression 
for this anomalous dimension in the case of n = 11. In principle, these calculations can be 
done with the help of the analytical methods of Ref.|[20||. 

In order to get rather rough estimate of the value of 7jvs,F 2 (l) we extrapolated the smooth 



curve, obtained by the interpolation procedures described above, to the point 



n 



1 and 



obtained 7^5_p 2 (l) = 306.7 for / = 4. Notice, that while the exact expressions for 7^ = 

and 7^5i? 2 (l) = 2.5576 do not depend from the number of flavors / taken into account, 

(2) 

the value 7^5. F2 (l) can be /-dependent. Indeed, repeating our interpolation/extrapolation 
procedure in the case of / = 3 we get with the non-controllable error-bars. Using this value 
as the input parameter in the renormalization-group expression for the first NS moment of F 2 
SF we arrive to the following order 0(a^) estimated expression for the Gottfried sum rule in 
the case of u — d flavor-symmetrical sea for / = 3 numbers of active flavors: 



GSR(Q 2 ) = 

where a s = a s (Q 2 
reads GSR{Q 2 ) = 



[F l 2 p { X) Q 2 )-F l 2 n {x,Q 2 ] 



, dx 



x 



1 + 0.036 



(I) 



0.72 



a. 



71 



(12) 



In the case of / = 4 numbers of flavours the analogous expression 



(1/3) 



Thus, taking into account the 



f 0.038(a s /7r) + 0.55(a s /vr) 2 

small positive QCD corrections to this sum rule, including the newly estimated order 0(a 2 )- 
contribution, can not explain the deviation of the theoretical prediction of Eq.(|12"D from the 
experimental result of the NMC collaboration, namely GSR = 0.235 ± 0.026 at Q 2 = 4 GeV 2 



E 



32| without the assumption, that the light quark sea is flavor asymmetric, which means that 



u - d > 0. 

4. The concrete fits of the available CCFR data of Ref.fll], made by means of the Jacobi 
polynomial expansion method, were performed by us for different values of the scale Qq with 
the help of two independently written computer programs. The fisrt Program has been created 
as the result of the works of Refs.|| [7j, while the second Program previously found its most 



distinguished application in the process of the fits of Ref. [23]. We have checked that both 



Programs are giving the identical results and thus are compatible with each other. 

In fact we confirmed the findings of Refs.[[| |7j that the results of the NS fits (namely 
the output values of and of the parton distributions A(Qq), b(Ql) ) c(Ql) 1 r y(Ql)) are non- 
sensitive to the choice of a, (3 provided large enough number of the Jacobi polynomial N max 
is taken into account (we are taking N max = 10). However, if smaller number of the Jacobi 
polynomial are considered, the reasonable convergence is achieved in the case of taking a ~ 0.7 
and 3 ~ 3. These values are allowing to approximate a posteriori with maximal precision the 
x-behavioral of the shape of the SF xF% by the weight function of the Jacobi polynomial in 
Eq.(ll) in the scaling limit and are determined by the values of the parameters b and c in the 
used model of the parton distributions (from the results of the previous Jacobi polynomial 
expansion fits of the CCFR data one can see that b(Q 2 = 10 GeV 2 ) = 0.717 ± 0.012 and 
c(Q 2 = 10 GeV 2 ) = 3.395 ± 0.043 @). We faced this problem of fixing sharply the values of 
a, 3 in the process of the NNLO fits, which we were not able to complete with the reasonable 
precision for N max > 6 in view of the uncertainties of the extrapolation for 7]vs( r at n = 
Nmax + 4 > 10. In principle, the question of the extension of our NNLO analysis to the case 
of N max = 7 can be studied in future if the exact expressions for 7^5^(11) become available. 

The results of the NLO fits are summarized in Table 4, where they are compared with the 
ones, obtained by the CCFR collaboration (see Ref.|L|). We have checked, that the results 
of our fits of the CCFR data are almost nonsensitive to the variation of Qq in the region 
3 GeV 2 < Ql < 50 GeV 2 . For definiteness, the results of Table 4 are presented in the case 
of taking Qq = 10 GeV 2 . Notice, that in spite of the fact that the methods applied in the 
process of our fits are different from the ones, used by the CCFR collaboration, the agreement 
with the results of Ref. [|]] is rather good. This feature supports further concrete applications 
of the Jacobi polynomial expansion method. 

Another important observation comes from the comparison of the values for A^-, presented 
in Table 4, with the outcomes of the previous fits of the available CCFR data for xF% SF by 
means of the Jacobi polynomial expansion method in the approximation when the TMC were 
neglected |§. This comparison indicates, that the inclusion of the TMC into the analysis 
lowers the obtained NLO values of A^^ by over 25 MeV. This observation confirms the 
previous findings of the IHEP-ITEP collaboration |34| that it is extremely important to take 
into account the effects of the TMC in the analysis of the neutrino-nucleon DIS data. 

The inclusion of the NNLO QCD corrections into the analysis of the available CCFR 
data for xF 3 SF made by us following the discussions presented above, result in the less 
vivid effect. In Table 5 we present the results of these fits for the value of Qq = 10 GeV 2 . 
One can see that for the maximal number of the Jacobi polynomials N max = 6, for which 
the reasonable convergence of Eq.(ll) was achieved, the resulting change (AA 



(A^) NNLO - (A^)iVLO ~ +2 MeV is very small. 

The results of the combined NS fits of xF 3 and F 2 data, which were obtained in the case of 
the cut Q 2 > 15 GeV 2 , considered in the works of the CCFR collaboration |], are presented 
at Table 6 in the case of taking Qq = 10 GeV 2 . It should be noted, that the combined 
analysis of xF 3 (x < 0.5) and F 2 (x > 0.5) CCFR data of Ref.|l| was made using the following 



MS 



I NNLO 



approximate conditions: first the TMC in Eqs.(9), (10) were taken into account separately to 
the NS moments M^f 3 c (Q 2 ) and M^ 2 C '(Q 2 ), but then both remaining moments M^ 3 (Q 2 ) 
and M^p 2 {Q 2 ) were considered to be equal to each other, parametrized at the scale Ql by 
the same way and evolved to the points Q 2 (where the experimental data is available) using 
the same renormalization-group equations with the same coefficient functions and anomalous 
dimensions as in the case of the renormalization group evolution of the Mellin moments of xF 3 
SFs (see Table 1 and Table 2). Then the Jacobi polynomial expansion of Eq.(8) was applied 
in the regions x < 0.5 and x > 0.5 with the same parameters a and (3. This approximate 
procedure has some theoretical support. Indeed, it can be shown, that in the region of large 
enough values of x, namely 0.5 < x < 1, the behavior of the reconstructed NS SF in Eq.(ll) 
is mainly determined by the Mellin moments M^ S (Q 2 ) with large enough n, for which the 
values of the factors A?(n)/4 = (C$(n) - C^(n))/4, A%\n)/16 = (C|?(n) - C^{n))/16 

and A 7 (n) = lNS,F 3 i n ) ~ lN ) s,F 2 ( n ) are ra t ner small. The applied approximation allows us not 
to double the number of the parameters of the fits, namely the Jacobi polynomials parameters 
a, (3 and the parameters of the parton distributions A, b, c and 7. 

However, in view of the violation of the approximations A^ (n) = 0, A^ (n) = 0, A 7 (n) = 
for the moderate values of n (see Table 1 and Table 2), this approximate procedure has 
definite theoretical uncertainties. Thus, to our point of view, the price for the decreasing 
the experimental uncertainties due to the inclusion in our fits of the F 2 CCFR data in the 
region of x > 0.5, which has higher statistical maintenance, is the increase of the theoretical 
uncertainties. The similar problem is also appearing in the course of the combined NLO fits 
of the CCFR data made in Ref. U using the approximation xF 3 = F2 at x > 0.5. Therefore, 
we think that in the summary tables of a s measurements it is better to present not only the 
results of the fits of the combined xF 3 and F 2 CCFR data, like it was done say in Ref. p5| , 
but the results of the separate fits of xF 3 CCFR data also. 

5. In order to study the structure of the perturbative QCD predictions for the NS 
moments of the SFs of the DIS neutrino-nucleon scattering, we performed the fits of the 
available CCFR data of Ref. [|1] at the LO also. The value for A^b turned out to be very stable 
to the variation of Q . For the cut Q 2 > 10 GeV 2 the result is = 131 ± 26 MeV in the 
case of the analysis of the xF 3 data, while in the case of the combined analysis of the xF 3 and 
F 2 data with the cut Q 2 > 15 GeV 2 we got A^b = 208 ± 26 MeV. Comparing these results 
with other ones presented in Tables 4-6 we arrive to the following conclusions 

• the effects of the NLO QCD corrections are important and are increasing the values of 
the QCD scale parameter by over 50 MeV in the case of 2^3 fit and by over 25 MeV 
in the case of the combined NS fits of xF 3 and F 2 data; 

• on the contrary to the findings revealed in the process of the NNLO NS analysis of 
the BCDMS data for F 2 SF of the muon-nucleon DIS, which was made in Ref. |[23|| , the 
inclusion of the NNLO QCD corrections into the analysis of the available CCFR data 
for xF 3 SF do not affect essentially the values of A-^=, extracted at the NLO level: 

• this effect is typical to the analysis of the data for xF 3 SF only. Indeed, after taking 
into account the CCFR data for F 2 we found that the inclusion of the NNLO QCD 
corrections into the analysis leads to the decrease of the value of A^- by over 25 MeV, 
like it was in the case of the NNLO NS fits of the BCDMS data for F 2 SF |3|; 

• we reveal the very interesting effect of the inclusion of the NNLO QCD corrections 
in the analysis of the CCFR data, namely the decrease of the difference between the 



17 



central values of the parameter A^L, which was extracted at the NLO from the separate 
NS fits of the xF 3 and xF 3 + F 2 CCFR data of Ref . [1] . Indeed, at the NLO the difference 
was over A£££ = (A^) xFs - (A^) xF3+F2 = -53 MeV, while at the NNLO it is shrunk 
by over 50% A and reads A^ NLO = -26 MeV (!!); 

J A MS V ' 

• another important effect of the inclusion of the NNLO corrections was revealed as the 
result of the fits of the available xF 3 CCFR data with the different Q 2 -cuts. We have 
found, that the difference between the values of A™-, obtained during the LO fits of the 
data with the cuts of Q 2 > 15 GeV 2 and Q 2 > 5 GeV 2 , which is over 28 MeV, is shrunk 
to 13 MeV in the NLO and to over 4 MeV at the NNLO (!!!). Thus we think, that 
in order to make the outcomes of the concrete fits less sensitive to the typical scales of 
the cuts of the experimental data it is very important to take into account NNLO QCD 
corrections. 

In order to extract from our analysis the related values of the QCD coupling constant a s (Mz) 
we first evolved the results presented in Table 5 and Table 6 through the threshold of the 
production of the 6-quark directly in the MS'-scheme using the derived in Ref.|]36|| relation 

/ \2/23 

Ag_ = AgU AgL/mJ Zri[(m b /Agy 2 ](- 963 / 13225 ) , where we chose m b « 4.6 GeV. After 

(5) 

this we substituted the numerical expressions for A^ into the NLO and NNLO inverse-log 
approximations for a s (taking / = 5 numbers of flavors and Q 2 = Mf ) and obtained the results 
presented in Table 7, where the systematical uncertainties were fixed using the estimates of 
the related errors in the published results of the CCFR collaboration |]] (see Eqs.(l),(2)). 
Looking carefully at the results of Table 7 we arrive to the following conclusions: 

• The previously omitted from the summary tables of Ref. [j35fl value of a s (Mz), extracted 
from the NLO analysis of the available xF 3 CCFR data, is in good agreement with 



the small value of a s (Mz) ~ 0.109, extracted in Ref. ||37|| from the QCD sum rules 
for the production cross section of bottomium states in e + e~-annihilation in the NLO 
of perturbative QCD. However, on the contrary to our analysis, the result of Ref.pTfl 
does not include the estimates of the possible effects of the NNLO perturbative QCD 
corrections. We think, that these effects can be important. 

Indeed, we observed, that the inclusion of the NNLO corrections into game makes the 
results of the fits of the available CCFR data more stable both to the procedure of 
neglecting or taking into account the data for F 2 SF at large x and to the typical scales 
of the cuts of the low Q 2 xF 3 data. 

Our final NNLO results are very closed to the result a s (Mz) = 0.113 ± 0.003(exp) ± 
0.00A(theor) [0q| , which comes from the NLO analysis of the BCDMS data by means of 
the integro-differential approach. In view of the appearance of the results of the NNLO 



calculations for the anomalous dimensions of the F 2 SF in the singlet case |17] it will be 



of interest to refine the NS NNLO Jacobi polynomial expansion analysis of the BCDMS 



data [23] in order to get the NNLO variant of the result of Ref.pq 



6. The discussions of the important questions of the theoretical uncertainties are now in 
order. 

It should be mentioned, first, that instead of using the inverse log expansion for a s (see 
Eq.(6)) one can solve the renormalization group equation for the QCD /5-function explicitly and 



arrive to the corresponding transcendental equation, used in Ref.|23[. The application of this 
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equation gives NLO analogs of the values of A-jiL, presented in Tables 5 and 6, systematically 
over 10 MeV smaller. However, we have checked that this difference is compensated by 
the differences, which appear in the process of the determination of the values of a s using 
the inverse logarithmic expression and the solution of the presented above transcendental 
equation. Thus, the final values of a s {Mz) turn out to be almost identical. 

The more important uncertainty comes from the arbitrariness in the choice of the matching 
point in the relation of Ref. [36]. (for the recent discussion of this problem see Ref . ] 39|) . It 
was mentioned in the detailed reviews of Ref.|30| that by varying the value of this point in the 
region Q m (0.75 — 2.5)m& one can estimate roughly an error of taking into account threshold 
effects. This estimate gives Aa s (M z ) = ±0.0015 PD| . Quite recently another interesting 
procedure of taking into account threshold effects was proposed |j41| . One can hope, that 
the continuing at present studies of the possibility of the application of this procedure in the 
NLO analysis of the CCFR data of Ref.0 for the xF 3 SF (see Ref.(42|) can give the more 
quantitative estimate of this theoretical uncertainty. 

Another interesting problem is related to the subject of taking into account the nuclear 
effects in the analysis of the data for the DIS of neutrinos on the targets with the fixed nuclear 
content. We should warn, however, that the NLO analysis of Ref. [O demonstrated that the 
effect of the influence of the nuclear corrections to the procedure of the extraction of the NLO 
value of A^ from the CCFR data of Ref. § for the xF( e SF can not be large and is decreasing 

A^ by < 7 MeV @. 

As follows from the discussions presented above, the effects of the NNLO perturbative 
QCD corrections are more vivid. Apart of the decreasing the difference between the values of 
a s (M z ), extracted from the NLO NS fits of xF 3 and xF 3 + F 2 data of Ref. |L|, the procedure 
of taking into account NNLO perturbative QCD corrections should decrease the theoretical 
uncertainty due to the arbitrariness in the choice of the renormalization scheme. This effect 
manifested itself most obviously in the process of the first detailed analysis of the previous 
CCFR data for the Gross-Llewellyn Smith sum rule ||44|| , which was made in Ref. |[45|| . The 



effect of the reduction of the scheme-dependence uncertainty at the NNLO was also confirmed 
in the process of the NS analysis of the BCDMS data |[23]. We think, that in our NNLO results 
of Table 7 this uncertainty can change the value of a s (Mz) by over Aa s (Mz) ~ ±0.0005, like 
it was in the case of the considerations of Ref. |I5|] . 

The value of this uncertainty can be simulated in part by taking into account the estimates 
of the non-calculable iV 3 LO-terms using the scheme-invariant methods, developed in the pro- 
cess of the related QCD studies for the DIS sum rules ffBfl. It is known from the considerations 



of Ref. [47| , that the results of applications of these methods are in agreement with the analo- 
gous results, provided by the Pade-resummation technique, which was previously used in the 
process of the fits of the DIS data of the BCDMS collaboration in Ref. ||4S|| . Therefore, it might 
be of interest to generalize the studies of the work of Ref. |^S[ to the case of the NNLO analysis 
of the CCFR data and thus to put the proposed above estimate of the scheme-dependence 
uncertainty of the NNLO value of a s (Mz) on the more solid ground. 

There are also some other yet non-fixed in our studies theoretical uncertainties, which 
come from the non- considered by us effects of the higher-twist corrections (which, however, in 
the considered kinematical region are expected to be less important than the included in our 
analysis effects of the TMCQ) and the effects due to non- zero masses of the c- and b-quarks. In 
view of this we prefer to follow the point of view of Ref. HO] that it is better to be more careful 
(than overoptimistic) in the process of estimating the total theoretical uncertainties of the 

5 This statement finds its empirical support in the results of the NLO QCD analysis of the neutrino- nucleon 
DIS data of the CDHS and BEBC-Gargamelle collaborations, which was made in Ref. [M some time ago. 
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values of a s (Mz)- Thus we are presenting the following NNLO results for a s (Mz), extracted 
from the available CCFR data of Ref. |l| by means of the Jacobi polynomial technique: 

xF 3 : a s {M z ) = 0.109 ± 0.003 (stat) ± 0.005{syst) ± 0.003 (theor) , (13) 

xF 3 + F 2 : a s (M z ) = 0.111 ±0.002 (stat) ±0.003 (syst) ±0.003 (theor) . (14) 

The presented theoretical errors are including all discussed above uncertainties. We are also 
aware, that the members of the CCFR group are planning to update their data of Ref.|], 
treating better the radiative EW corrections, dimuon and charmed mass corrections, along 
with extending kinematical coverage and reducing systematical errors |5(J . These studies can 
affect the values of the available CCFR data at low x region (the necessity of this improvement 
was previously emphasized in Ref.|^T|, devoted to the extraction of the parameters of the 
CTEQ and MRS parton distributions). In spite of the fact that in our analysis we have cutten 
down the low Q 2 and thus low x points, it can be of real interest to perform the more detailed 
NNLO analysis of the updated CCFR data in future, after they will become available. 

Another interesting problem, which is more closely related to the extraction of the value 
of the Gross-Llewellyn Smith sum rule from the extrapolated data of the CCFR collabora- 
tion H |3B|, is the determination of the values of the parameters of the parton distributions 
A(Ql),b(Ql),c(Ql)^(Ql) for the SF xF 3 from the LO, NLO and NNLO fits. In view of the 
lack of space, we are postponing the presentation of the concrete results of our fits up to 
another publication. Note, however, that in spite of the existence of the results of the recent 
detailed studies of Ref. |52| we are still doubting whether this simple parametrization can be 



applied with the necessary theoretical precision in the region of very small x. We hope, that 
the possible future theoretical considerations of the description of the behavior of the NS 
structure functions at very small x, which might be done, say, using the developed in Ref. ||53|| 
methods, can be useful to clarify this interesting problem. 

7. In conclusion, we would like to stress again, that using the Jacobi polynomial expansion 
method, developed in Refs.|5], || [?|, we have shown that the effects of the NNLO corrections 
are important in the analysis of the most precise up to now experimental data of the CCFR 
collaboration for the NS SFs of the neutrino- nucleon DIS. The obtained NNLO values of 
a s (Mz), extracted separately from the fit of the data for the xF 3 SF and from the combined 
NS fit of xF 3 + F 2 data, turn out to be in good agreement with the previous NLO determination 
of a s (M z ) from the fit of the BCDMS data for F 2 SF of the muon-nucleon DIS |3§. Thus the 
problem of understanding of the existing differences between the values of a s (Mz) extracted 
from DIS and from the combined fits of the data of 4 LEP groups for the Z° — > hadrons 
decay width is still on the agenda. For the certain attempts to understand this discrepancy 
by introducing the effects of virtual SUSY particles and other effects of the possible new 
physics (see e.g. Ref.p>4| and Ref.[^5|). Note, however, that there are certain systematical 



discrepancies between the values of a s (Mz), extracted from the Z° — > hadrons decay width 
by the DELPHI, ALEPH, OPAL and L3 groups separately (see e.g. [f)6|]). Therefore, it might 
be also of interest to understand better the origin of the possible existing discrepancies. 
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Table 1. The numerical values of the NLO and NNLO coefficients for the coefficient functions 
of the NS moments of F 2 and xF 3 SFs at / = 4 number of flavours. 



n 
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2.5576 





306.6810 


2 


71.3745 


71.2410 


612.0598 
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100.8013 


100.7819 


837.4264 
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120.1447 


120.1401 


1005.8235 
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134.9049 


134.9035 


1135.9189 
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147.0029 


147.0024 


1242.0056 
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157.3323 


157.3321 


1334.0017 
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166.3862 


166.3861 


1417.4506 
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174.4683 


174.4682 


1493.5205 
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181.7808 


181.7808 


1559.0048 
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188.4662 


188.4662 
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194.6293 


194.6293 
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13 


200.3496 


200.3496 
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14 


205.6891 


205.6891 


? 



Table 2. The used numerical expressions for the NLO and NNLO coefficients of the anomalous 
dimensions of the moments of the NS SFs at / = 4 number of flavours. 
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lNS,F 2 (n) 


A 7 


A 7( n )/7S,F 2 H 


2 


71.3745 


1.335 x 10" 1 


1.87 x 10" 3 


4 


120.1447 


4.6 x 10" 3 


3.83 x 10" 5 


6 


147.0029 


5 x 10" 4 


3.4 x 10~ 6 


8 


166.3862 


1 x 10" 4 


6.01 x 10~ 7 



Table 3. The related numerical differences between the NLO coefficients of the anomalous 
dimensions of the NS moments of F 2 and xF 3 SFs at / = 4 number of flavours. 



1 c 



Q 2 cut (GeV 2 ) 


x cut 


Order 


A^L (MeV) 

MS V 1V1C v J 


y 2 ( xFi) /neD 


O 2 > 5 


r < 7 


NLO 


1 78 ± 29 


81 6/81 




CCFR result 




170 ± 31 


83 8/81 

. w / OX 


O 2 > 10 


r < 7 


NLO 


182 ± 30 

X UZj _l "XV/ 


64 5 /65 




CCFR result 




171 ±32 


66.4/65 


Q 2 > 15 


x < 0.7 


NLO 


191 ±34 


51.4/54 




CCFR result 




179 ±36 


53.5/54 


Q 2 > 15 


F 2 x > 0.5 


NLO 


235 ± 28 


51.9/54 




CCFR result 




210 ±28 





Table 4. The comparison of the results of our NLO fits of xF 3 CCFR data of Ref. [1] with 
the results provided by the CCFR collaboration in Ref. [1] (x 2 is normalized to the number of 
experimental points (nep).) 



Q 2 cut (GeV 2 ) 


Order 


N 

1 y max 


(MeV) 


X 2 (a;F 3 )/nep 


Q 2 > 10 


LO 


10 


131 ±26 


76.7/65 


Q 2 > 10 


NLO 


10 


182 ±30 


64.5/65 


Q 2 > 10 


NNLO 


6 


184 ±31 


57.6/65 



Table 5. The results of our NNLO analysis of the CCFR data of Ref.[l] for xF 3 SF. 



Q 2 cut (GeV 2 ) 


Order 


N 


(MeV) 


X 2 (xF 3 )/nep 


Q 2 > 15 


LO 


10 


208 ± 26 


62.3/54 


Q 2 > 15 


NLO 


10 


235 ± 28 


51.9/54 


Q 2 > 15 


NNLO 


6 


210 ±24 


44.2/54 



Table 6. The results of our combined NNLO NS analysis of the CCFR data of Ref.[l] for 
xF 3 and F 2 SFs. 



Order 




a s (M z ) 


Aa stat 




NLO 


xF 3 


0.108 


0.003 


0.005 


NLO 


xF 3 + F 2 


0.112 


0.002 


0.003 


NNLO 


xF 3 


0.109 


0.003 


0.005 


NNLO 


xF 3 + F 2 


0.111 


0.002 


0.003 



Table 7. The obtained NLO and NNLO results for a s (M z ). 



1 c 



